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ABSTRACT 

We investigate the non-thermal X-ray emission from Cas A, using BeppoSAX, CGRO-OSSE, and 
Chandra data. For the hard X-ray continuum we test the model proposed by Laming, which invokes non- 
thermal bremsstrahlung from electrons accelerated by lower hybrid plasma waves. The justification for 
this model comes from our determination of a lower limit to the average magnetic field oi B > 0.5 mG. For 
such high magnetic fields the synchrotron losses are severe enough that most of the electron populations 
responsible for the radio emission have maximum electron energies well below the limit for which X-ray 
synchrotron emission is important. However, we do suggest that the rim surrounding Cas A, seen in 
Chandra continuum images, is X-ray synchrotron emission. The width of this rim of 1.5" to 4", can 
be used to infer the magnetic field near the shock front, for which we estimate B = 0.08 — 0.16 mG, 
and electron energies of 57 — 40 TeV. This magnetic field strength is lower than the average magnetic 
field, but higher than what may be expected from shocked interstellar medium, suggesting either a high 
magnetic field in the wind of the progenitor, or rapid, post shock, field amplification by non-linear growth 
of plasma waves. Combining the two magnetic field measurements we have constructed a simple two 
zone model. Most of the radio emission comes from inside Cas A where the magnetic field is strong. In 
contrast, the inverse Compton emission is dominated by emission from near the shock front. Only for 
our lower limit on the magnetic field strength near the shock front is it possible to explain the recent 
detection of TeV emission by HEGRA with inverse Compton, for which, in addition, we have to assume 
a rather high far infrared photon density that should be twice as high as our best estimate of ^ 70 cm""^. 
Pion decay is therefore likely to be the dominant emission from Cas A at TeV energies. 

Subject headings: shock waves ~ magnetic fields - gamma rays: observations - X-rays: observations 
individual (Cassiopeia A) - supernova remnants - 



1. INTRODUCTION 

Supernova remnants have long been thought to be the 
source of cosmic rays with energies up to at least 10"'^^ eV, 
where there is a break in the observed cosmic ray spec- 
trum, usually referred to as the "knee" . The recent dis- 
coveries of X-ray synchrotron and TeV emission from some 
supernova remnants support this long held belief by con- 
firming that electrons are accelerated to energies of at least 
a few TeV (Reynolds & Keohane 1999). This is still some- 
what short of the "knee" , but unlike for ions, the maximum 
electron energy may be limited by radiation losses. 

SN 1006 has become the canonical example of a su- 
pernova remnant with X-ray emission dominated by 
synchrotron radiation from shock accelerated electrons 
(Koyama et al. 1995). The presence of electrons in excess 
of 10 TeV is further supported by the detection of TeV 
emission with the CANGEROO Cerenkov detector (Tani- 
mori et al. 1998). The observed TeV emission is most likely 
inverse Compton scattering by the same population of elec- 
trons that is responsible for the X-ray continuum (Tani- 
mori et al. 1998), although there is some debate about this 
interpretation (Berezhko et al. 2002). 

Another often mentioned candidate for an X-ray syn- 

^ Chandra fellow 



chrotron emitting supernova remnant is Cas A. Its soft 
X-ray emission is not continuum dominated, but the ob- 
served emission above 20 keV cannot be explained by ther- 
mal bremsstrahlung, and has therefore been attributed to 
synchrotron radiation (The et al. 1996; Allen et al. 1997; 
Favata et al. 1997). The detection of TeV emission by 
HEGRA is consistent with this interpretation, but pion 
decay associated with relativistic ions is an equally valid 
interpretation (Aharonian et al. 2001). 

For Cas A the interpretation of both the hard X-ray 
and the TeV emission is critically dependent on the mag- 
netic field strength. Equipartition arguments to explain 
the radio luminosity imply magnetic fields oi B ^ 0.4 - 2 
mG (Longair 1994; Wright et al. 1999), far in excess of 
the canonical value for the interstellar medium of 3 /xG. 
The typical loss time of an electron with energy E, due to 
synchrotron radiation is (e.g. Reynolds 1998): 



For a magnetic field of 1 mG, a particle with E — 8 TeV 
loses its energy in 2 yr. For Cas A, 8 TeV is the max- 
imum energy needed to explain the hard X-ray emission 
(Reynolds & Keohane 1999). The two year lifetime is much 
shorter than the lifetime of Cas A, which is probably the 
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remnant of a supernova that exploded in AD 1680 (Ash- 
worth 1980). Similarly a high magnetic field lowers the 
number density of relativistic electrons for a given radio 
brightness, which in turn implies a low inverse Compton 
contribution to the observed TcV omission. 

One way out of this situation would be to assume that 
electrons are continuously being reaccelerated, but that 
seems unlikely. The radio morphology of Gas A shows a 
bright ring of radius ~1.8', whereas the shock front, the 
primary site for cosmic ray acceleration, is at a radius of 
-^2.6'. Moreover, detailed radio spectral index studies do 
not find evidence of reacceleration in bright radio knots, 
which arc probably cjccta fragment penetrating the shell 
(Anderson & Rudnick 1996; Wright et al. 1999). These 
studies rather suggest that brightness variations are the 
result of variations in the magnetic field, which light up 
the diffuse relativistic electron population in various de- 
grees. 

Here we test an alternative interpretation of the 
hard X-ray emission with BeppoSAX-PDS data, namely 
bremsstrahlung from electrons accelerated by lower hy- 
brid waves, which was suggested by Laming (2001a, b). 
An additional argument for testing such a model is that 
XMM-Newton observations show that even above 10 keV 
a substantial fraction of the continuum emission is not as- 
sociated with the forward shock (Bleeker et al. 2001). The 
BeppoSAX data, in combination with CGi^O-OSSE data, 
also provides an upper limit on the bremsstrahlung from 
the electron cosmic ray injection spectrum, from which a 
new lower limit on the average magnetic field strength can 
be obtained. 

Nevertheless, we argue that X-ray synchrotron radiation 
does exist in Gas A, and is coming from close to the for- 
ward shock, where it shows up as a narrow rim surrounding 
the remnant, observed in Chandra X-ray continuum im- 
ages (Gotthelf et al. 2001). Interpreting the width of the 
rim as a typical length scale for synchrotron losses, we are 
able to constrain the magnetic field near the shock front. 
Our analysis supports the idea that the average magnetic 
field is close to its equipartition value, whereas we find 
that near the shock front the magnetic field is lower, but 
still an order of magnitude higher than what may be ex- 
pected from shocked interstellar medium (~ 10~^G). We 
will discuss the implications of these values for the hard 
X-ray and TeV emission. 

2. DATA AND OBSERVATIONS 

Our analysis of the X-ray emission of Gas A is based 
on BeppoSAX-Vm data and archival CGi?0-OSSE and 
Chandra data. The BeppoSAX-FT)S (Frontera et al. 1997) 
and CGi20-0SSE (Johnson et al. 1993) instruments both 
consist of four phoswitch detectors (Nal(Tl) /GsI(Na) scin- 
tillation detectors), behind two rocking collimators. They 
differ, however, substantially in energy range, field of view, 
and effective area. The BeppoSAX-VDS instruments cover 
an energy range of 15-300 keV, whereas this is 50 keV to 
10 MeV for CGi?0-OSSE. 

BeppoSAX observed Gas A for 500 ks in May and June 
2001 in order to observe the 67.9 keV and 78.4 keV nuclear 
decay lines of ^''Sc, the nuclear decay product of ''^Ti. The 
detection of these lines at the 3-5c7 level, depending on the 

^ In addition to the observations listed in The et al. (1996), this me; 



hard X-ray continuum model, was reported by Vink et al. 
(2001), which provides more details about the data set. 
Due to the rocking collimators the effective exposure time 
is about half the observation time. Together with archival 
BeppoSAX data the spectra used for our analysis here is 
based on a total effective PDS exposure of 311 ks. 

The spectrum, as presented in Vink et al. (2001), clearly 
shows the hard X-ray spectrum with a bump associated 
with the ^^Sc line emission. Marginal {2a) residual emis- 
sion around 60 keV is possibly caused by fluorescent tan- 
talum K shell emission, as a result of cosmic rays interact- 
ing with the collimator, which consists of tantalum with 
a thin copper-tin bi-layer. Some residual emission above 
100 keV is not easily explained (cf. Fig. 1 & 2). It may be 
intrinsic to Gas A, but some instrumental contamination 
seems more likely. The residual emission peaks around 
200 keV, but is not narrow enough to be caused by fluo- 
rescent line emission from any particular material. There 
are no other similar deep BeppoSAX observations with a 
moderate hard X-ray brightness like Gas A, which would 
allow us to identify possible instrumental features at this 
level. 

Although the residual emission above 100 keV is statis- 
tically not very signiflcant, one of our concerns is to obtain 
an upper limit on the emission above 100 keV in order to 
obtain a lower limit to the average magnetic field in Gas A. 
We therefore included archival GGi^O-OSSE spectra in our 
analysis. Moreover, the effective area of the instruments 
peaks around 250 keV, making it useful for constraining 
the emission of Gas A above 100 keV. Gas A was observed 
several time by GGi?0-OSSE (see The et al. 1996). We 
added all available Gas A spectra and averaged the corre- 
sponding response matrices. The observations were made 
during the viewing periods 34 - 815 with a total exposure 
of 909 ks.^ The total exposure time was 909 ks. 

In addition to hard X-ray data, we use a deep (50 ks) 
archival Chandra observation with the AGIS-S3 chip to 
complement our analysis. The observation was made in 
January 2000. The same data were used by Hwang et al. 
(2000). For our analysis we used level 2 archival data prod- 
ucts, generating spectra and calibration products with the 
standard Chandra software package Ciao vl.1.0. 

3. THE LOWER HYBRID ELECTRON ACCELERATION 

MODEL 

Gas A's morphology indicates a complicated, turbulent 
hydro dynamical structure. Radio, optical and X-ray im- 
ages all show many bright knots on top of more diffuse 
emission. As the forward or reverse shocks strike and 
illuminate these features, secondary shocks will split off 
from these primary shocks and propagate back into the 
already shocked circumstellar or ejecta plasma. In this 
way shocks may propagate across the likely region of high 
magnetic field at the contact discontinuity, where field am- 
plification by Rayleigh- Taylor instability may occur. The 
acceleration of electrons by lower hybrid waves generated 
by these shocks has been considered by Laming (2001a, b). 
Bremsstrahlung emitted by these accelerated electrons was 
shown to provide a very good match to the hard X-ray con- 
tinuum of Gas A, using the ^ 50 ksec of BeppoSAX-PT)S 
data that were available at the time. 

VP 617, 715, 810, 814, and 815. 
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Lower hybrid waves (LHW) are electrostatic ion oscil- 
lations directed within an angular range zbwpi/ojpe of the 
direction perpendicular to the magnetic field, where 
and ujpi are the electron and ion plasma frequencies respec- 
tively. The electron screening that would usually damp 
such oscillations is inhibited by the magnetic field, giving 
a criterion on the electron temperature that the electron 
gyroradius be smaller than the wavelength divided by 27r. 
Since uj/k± << a;/fc|| the wave can simultaneously be in 
resonance with ions moving across the magnetic field lines 
and electrons moving along them. Consequently energy 
equilibration between ions and electrons may proceed at 
a much faster rate than by Coulomb collisions alone, and 
saturate when the electron thermal pressure is equal to 
the magnetic pressure. Such processes have been consid- 
ered in the context of two-temperature accretion flows by 
Begelman & Chiueh (1988). Here cross field ion motions 
are produced by a combination of curvature and diamag- 
netic drift, and the magnetic field is presumably gener- 
ated by a magneto-rotational instability (Balbus & Haw- 
ley f99f, f998; Quataert et al. 2002), the critical issue 
being whether sufhcient electron-ion equilibration can oc- 
cur to destroy ADAF or similar solutions for the accretion 
flow. In dilute plasmas where the electron acceleration 
time is much shorter than the electron-electron collision 
time, a non-Maxwellian electron distribution in the com- 
ponent of velocity along the magnetic field can result. This 
has been modeled numerically using particle-in-cell simu- 
lations by Shapiro et al. (f 999) for application to cometary 
X-rays, and analytically by Vaisberg et al. (1983) and Kra- 
nosel'skikh et al. (1985). Lower hybrid waves generated by 
a modified two-stream instability and the resulting acceler- 
ated electrons were all observed at Comet Halley (Klimov 
et al. 1986; Gringauz 1985). Bingham et al. (2000) consid- 
ered LHW in supernova remnants and Laming (2001a, b) 
applied the analytic formulations, which agree well with 
the numerical work, to the case of accelerated electrons in 
Cas A. The non-Maxwellian electron distribution function 
depends on two parameters, the overall normalization (i.e. 
the fraction of plasma electrons that are accelerated) and 
the maximum electron energy in the distribution. In our 
fitting we treat these both as free parameters. A demon- 
stration of how these variables are related to other physics 
of the Cas A plasma, and hence could be determined is 
given in the Appendix. 

We fitted the predicted spectrum, calculated as in Lam- 
ing (2001a), but with a relativistic bremsstrahlung cross 
section taken from Nozawa, Itoh, & Kohyama (1998), to 
the observed BeppoSAX-PDS X-ray spectrum. Apart from 
the non-thermal bremsstrahlung spectrum several other 
components, which are likely to contribute to the hard X- 
ray emission, are included in the spectral model. At the 
low end of the hard X-ray spectrum, the tail of the the 
thermal bremsstrahlung with Tc = 4 x 10^ K is expected 
to contribute to the spectrum. In addition, ^^Sc line emis- 
sion at 68 and 78 keV is expected, as a result of the decay 
of '''^Ti synthesized during the explosion (Vink et al. 2001; 
lyudin et al. 1994). The line fluxes of the two nuclear tran- 
sitions are expected to be almost identical. A broadening 
corresponding with a velocity dispersion of 5000 km s~^ 
was included. A residual feature around 60 keV in the ob- 

See the X-ray data booklet (Thompson et al. 2001). 



served spectrum may be attributed to cosmic ray induced 
fluorescence from the tantalum collimator, and was mod- 
eled by delta functions corresponding to Ta Ka (57 keV) 
and K/3 (65 keV) with a ratio of 4:1."^ The spectrum and 
best fit model are shown in Fig. 1 and model parameters 
are given in Table Al. 
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Fig. 1. — BeppoSAX-PDS spectrum of Cas A with the best fit 
*'*Sc/non-thermal bremsstrahlung model. The individual emission 
components are: ''''Sc line emission (red), non-thermal continuum 
(blue), thermal continuum (green) and possible line contamination 
from collimator material (tantalum, cyan). The observed count rate 
in each channel has been divided by the effective area in order to 
yield approximate flux densities. 

The parameter values indicate that electrons are possi- 
bly accelerated up to ~95 keV. The model gives a very 
good fit to the spectrum up to 100 keV. However, it does 
not describe the emission observed above 100 keV. As dis- 
cussed above, the data above 100 keV may be less reliable, 
and the low flux point around 100 keV suggest that the 
actual continuum emission is indeed curved, as predicted 
by the LHW model. 

In Table Al we also list for comparison the parameters 
for a simple power law continuum spectrum, and the X-ray 
synchrotron models available in the spectral fitting pack- 
age xspec (Arnaud 1996), srcut and sresc (Reynolds 1998; 
Reynolds & Keohane 1999). The srcut model is a very 
simplified model of X-ray synchrotron emission from a rel- 
ativistic electron distribution with an exponential cut-off. 
In reality the cut-off energy is likely to vary from place 
to place, giving a more gradual decrease in slope. More- 
over, as explained in the introduction, for the bulk of the 
relativistic electrons, at a remote distance from the shock 
front, the maximum electron energy is probably far below 
what is needed to produce X-ray synchrotron emission. 
The sresc model takes these spatial effects into account, 
assuming a Sedov model for the evolution of the supernova 
remnant. Although, sresc is more physical, it may not be 
readily applicable to Cas A, which has large magnetic field 
gradients and is not well described by the standard Sedov 
model (Vink et al. 1998). 

The LHW model is in a similar way somewhat oversim- 
plified, as also here it is likely that maximum energies vary 
within the remnant, depending on local conditions such as 
the local shock velocity, obliquity and ion composition. 
Nevertheless, given the limited number of parameters, the 
model fits very well. Note that the spectral slope between 
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20-50 keV, is a property of the model, and not a free pa- 
rameter. 

Finally, note that using the LHW model results in 50% 
higher ^"^Sc line fluxes, than using a power law contin- 
uum, a result of the continuum curvature predicted by 
both models. The higher ^^Sc line flux is close to the lat- 
est CGRO-COMPTEL measurements of the '^^Cd. line flux 
(Schonfelder et al. 2000), and translates into a synthesized 
'''^Ti mass of 1.8x10"* M©, for a distance of 3.4kpc (Reed 
et al. 1995) and a *'*Ti decay time of 85.4 yr (Ahmad et 
al. 1998; Gorres et al. 1998; Norman et al. 1998). 

4. A NEW LOWER LIMIT TO THE AVERAGE MAGNETIC 

FIELD 

The LHW model is attractive to consider for Cas A be- 
cause of its high magnetic field strength, which is esti- 
mated to be '--'1 mG, based on the assumption of equipar- 
tition of the relativistic particle population and the mag- 
netic field energy (e.g. Longair 1994). Although there is 
no compelling physical reason why the equipartition as- 
sumption should be valid, estimates based on the absence 
of gamma-ray emission from Cas A supports a value of 
the magnetic field in excess of what might be expected 
from a shock compression of the interstellar magnetic field 
(Cowsik &; Sarkar 1980). The most recent estimate of the 
magnetic field is based on the upper limit obtained by 
CGRO-EGRET, suggesting B > 3.5 x IQ-^ G (Esposito 
et al. 1996; Atoyan et al. 2000). This value is derived 
by comparing the radio synchrotron flux density with the 
expected bremsstrahlung from the relativistic electrons re- 
sponsible for the radio emission. 

Here we use essentially the same method, but less di- 
rectly, as we use an extrapolation of the relativistic elec- 
tron population to energies observed by BeppoSAX-PDS 
and CGRO-OSSE. The extrapolation used is similar to the 
one by Asvarov et al. (1990), and assumes that the elec- 
tron momentum spectrum is a power law distribution, as 
is expected for an electron population generated by flrst 
order Fermi acceleration (Bell 1978). This means that the 
energy spectrum flattens at non-relativistic energies. The 
electron cosmic ray injection spectrum is not well known, 
but models indicate that at lower energies the spectrum is 
steeper than expected from extrapolating the relativistic 
spectrum (Bykov & Uvarov 1999). In that case our upper 
limit on the cosmic ray normalization is too high, which 
translates into too low a lower limit. 

Of more concern is the role of energy losses on the shape 
of the spectrum, as it will tend to flatten the low energy 
end of the spectrum. We checked the signature of adia- 
batic losses on the spectrum, and found that it does affect 
the spectrum mostly at mildly relativistic energies, where 
it results in a steeper spectrum at low energies, is flat- 
ter at intermediate energies and then assumes the original 
power law index at extreme relativistic energies. We did 
not take the steepening into account, as we do not know 
how much the average electron spectrum is affected by it. 
But if there is a steepening, it makes our estimates too con- 
servative. The most dominant other loss factor between 
100 and 1000 keV in the ionized plasma of the remnant is 
electron-electron collisions. For an electron temperature 
of 4 X 10^ K and density of 20 cm~^, the loss time of a 
100 keV electron is ~200 yr, rapidly increasing to ~500 



yr for 200 keV (Huba 1994). These loss times are compa- 
rable to or larger than the average plasma age of Cas A, 
indicating that non-adiabatic losses are likely to be small. 

The synchrotron flux density for a relativistic electron 
power law distribution, with energy index q, is given by 
(Blumenthal & Gould 1970): 



= 1.70 X 10 



21 / N ^'+1 ^6.26 10^^,-1 



where k is the particle normalization, V the emitting vol- 
ume, d the distance, B the magnetic fleld, and a{q) is a 
factor tabulated in Blumenthal & Gould (1970). 

Bremsstrahlung from the same population of electrons, 
which flattens at non-relativistic energies, scales with the 
factor T,iniZfKV/4TT(P, where TiiUiZ"^ indicates the sum of 
densities of all the ion species with charge Zi. So from the 
combination of synchrotron emission and bremsstrahlung, 
and an estimate of J^iUiZf, we can estimate k and B. 

The radio spectrum of Cas A has a spectral index of 
a = —0.78 and based on a compilation of radio data (e.g. 
Baars et al. 1977) we estimate a radio flux density of 2522 
Jy at 1 GHz for the epoch 2000. The spectral index in- 
dicates an electron energy spectrum with q ~ 2.56 at rel- 
ativistic energies and flattening to 1.77 at non-relativistic 
energies, resulting in an X-ray photon index of -2.2. This 
is at odds with the observed spectral index of -3.3, which 
means that 100 keV must have a different origin, taken 
here to be bremsstrahlung from LHW accelerated elec- 
trons. We adopt the LHW model, and try to obtain an 
upper limit on the emission from an additional compo- 
nent with a power law slope of -2.2 from the combined 
BeppoSAX und CGRO-OSSE data (Fig. 2). 
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Fig. 2.— Best fit model that combines the LHW/^^Sc with the 
electron cosmic ray injection spectrum. Shown are the BeppoSAX- 
PDS (black) and CGRO-OSSE (red) spectrum of Cas A. 

The best flt normalization for such a power law spec- 
trum is (4.0 ± 1.5) X lO"'^ ph s-^keV-^cm-^ at 100 keV, 
with a 2a upper limit of 6.2 x 10~^ phs~^keV"^cm"^. 
We measured the electron spectrum normalization more 
directly by fitting a model of bremsstrahlung from a rela- 
tivistically correct power law momentum spectrum, using 
the gaunt factors proposed by Haug (1997). The upper 
limit on the bremsstrahlung normalization obtained this 
way is K'EiniZfV/AncP < 65. Soft X-ray observations im- 
ply J:^n,Zf ~ 20 (Vink et al. 1996; WiUingale et al. 2001b), 
which, combined with the bremsstrahlung normalization, 
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gives for the average magnetic field B > 0.5 mG.^ This is 
higher than the value derived by Atoyan et al. (2000) from 
the CGRO-EGRET upper limit, especially taking into ac- 
count that they assumed a higher value for T,iniZf. Note 
that the above lower limit is not a volume average, but 
based on the ratio < > / < nY.imZ'^ >. The 

derived lower limit is in good agreement with equipartition 
of magnetic and cosmic ray energy. 

5. THE MAGNETIC FIELD NEAR THE SHOCK FRONT 

The high average magnetic field in Gas A is inconsistent 
with simple shock compression of the interstellar magnetic 
field, and indicates that the magnetic field inside the rem- 
nant must have been enhanced by turbulence. A likely 
place for turbulence is the contact discontinuity, which 
may explain why the radio emission has a bright ring 
of emission, as this is the region in the remnant where 
magnetic field amplification is important. A similar argu- 
ment was used by Anderson & Rudnick (1996) to explain 
why radio knots have the same radio spectral index as the 
surrounding diffuse emission; the enhanced magnetic field 
lights up the background electron population. One may 
wonder then, what the initial interstellar magnetic field 
was, before turbulent magnetic field amplification. 

An estimate can be derived from the observation of the 
shock front of Gas A by Chandra. In the X-ray contin- 
uum band from roughly 4.5-6 keV Gas A seems bound by 
a narrow rim of emission, which marks the onset of ra- 
dio emission. The narrowness of the rim indicates strong 
limb brightening of a layer, which much have some special 
emission properties. The fact that the radio and X-ray 
morphology are dissimilar as far as this rim is concerned, 
seems to argue against a synchrotron origin for the X- 
ray emission, see e.g. Gotthelf et al. (2001). But this is 
not true, if the width of the rim is a result of the limited 
lifetime of the ultra-relativistic electrons, as soon as they 
are swept away from the shock front. Further inside the 
shock the maximum energy of the electrons has so much 
decreased that no X-ray synchrotron emission is produced 
anymore. This is not the case for the radio emission, 
caused by electrons with much lower energies, for which 
synchrotron and inverse Gompton losses are negligible (e.g. 
Wright et al. 1999). This behavior can be seen in the mod- 
els by Reynolds (1998), in which the synchrotron emitting 
shell is thinner at higher frequencies. The steady increase 
in radio emission toward the bright ring of emission may 
be the effect of accumulation of relativistic electrons swept 
away from the shock front, and a strong gradient of the 
magnetic field. 

In this interpretation most of the electron populations, 
which reside outside the X-ray rims, have a maximum 
energjr cut-off well below 1 TeV, unless reacceleration is 
important. However, radio spectral index measurements 
do not indicate any reacceleration of electrons by ejecta 
knots, which would likely lead to different indices of radio 
knots and diffuse emission (Anderson & Rudnick 1996). 
Note that this argues also against the use for Gas A of the 
simple maximum curved synchrotron model employed by 
Reynolds & Keohane (1999), which merely assumes that 
all electron populations have the same energy cut-off, and 

This replaces the value reported in Vink & Laming (2002), in which 
Ginzburg & Syrovatski 1965). 



the magnetic field is constant. 
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Fk;. 4.- - Chandra ACIS-S3 X-ray spectra of the narrow rim 
of emission (lower spectrum) and from emission just inside the rim 
(xlO). Note the lower equivalent width of the line emission from the 
spectrum of the rim (cf. Table A2). 

The rim is much narrower than the ~l/12th shock ra- 
dius expected for plasma shocked with a factor 4 compres- 
sion ratio. That they are indeed consistent with X-ray 
synchrotron emission can be exemplified by a comparison 
of the X-ray continuum flux and radio flux of part of the 
rim in the northeast (see Fig. 3, and Table A2). We used 
archival Chandra data to obtain a spectrum from this fil- 
ament and determined that the X-ray continuum flux be- 
tween 4 and 6 keV is 1.0x10-^2 ergs-icm-2 (2.1x10"'^ Jy 
at 1.2 X 10^* Hz). We normalized an archival VLA radio 
map (epoch 1987) to a flux of 2522 Jy and extracted a 
flux density from the same region. As the shock front may 
have moved between 1987 and 2001 by as much as 5", 
we extracted also a flux from the same region, but moved 
5" toward the center of the remnant. This yielded a flux 
density estimate between 3 Jy and 6 Jy at 1 GHz, corre- 
sponding to a broad band spectral index of a = —1.15 (the 
factor two uncertainty in flux density is irrelevant over a 
frequency range of 9 decades). This is steeper than the 
radio spectral index of 0.78 and thus indicates that the 
X-ray flux density is below the extrapolation of the radio 
spectrum; consistent with the interpretation that part of 
the X-ray emission is synchrotron emission. 

Indeed, not all X-ray emission is synchrotron, as the nar- 
row rim displays line emission from Mg, Si and S, which 
should be accompanied by bremsstrahlung as well. The 
line emission from the rim is, however, weaker and the 
spectrum harder, than for the spectrum of a region just 
inside of the rim (Table A2). For example fitting a ther- 
mal non equilibrium model with an additional power law 
gives a spectral index of -2.2, whereas away from the rim 
the spectral index is unconstrained, but steeper than -4. 
Moreover, it is hard to find another interpretation for the 
rim than synchrotron emission, as only energy losses can 
explain the narrow features. For a hot plasma the losses 
are small and the gradual ionization and possible equili- 
bration of electrons and ions behind the shock would lead 
to a steady increase in intensity behind the shock, not in 
a sudden onset of emission, followed by a rapid decline. 

we missed a factor 47r in the synchrotron normalization (eq. 2, cf. 
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Fig. 3. — Detail of a VLA radio image (left) and Chandra-AClS X-ray image from a compilation of continuum dominated bands between 
2.7 - 9 keV. The spectral extraction regions are indicated in the Chandra image (NE rim and NE inside rim). 



Interestingly, identifying the rims with X-ray syn- 
chrotron emission from recently accelerated electrons 
makes this aspect of Cas A very similar to SN 1006, where 
the X-ray synchrotron emission is coming from a simi- 
larly narrow rim of emission, as shown by archival Chan- 
dra data. For SN 1006, however, the continuum emission 
from the rims are dominating the total continuum emis- 
sion, whereas for Cas A the X-ray continuum morphology 
shows that a substantial fraction of the X-ray continuum 
is associated with the bright ejecta shell. In SN 1006 the 
synchrotron emission comes mainly from the northeast- 
ern and southwestern limb, which is usually interpreted 
as caused by the large scale orientation of the magnetic 
field (Reynolds 1998). For Cas A, however, we do not see 
such axi-symmetry. This implies a much more randomly 
oriented magnetic field in the surrounding medium. 

Assuming that the limb brightened rim emits indeed 
X-ray synchrotron radiation, we can derive from the typ- 
ical width of the rim the magnetic field behind the shock 
front and the typical electron energy responsible for the 
synchrotron radiation. The electron energy may be asso- 
ciated with the typical exponential cut-off energy of the 
electron spectrum, but the very nature of the losses make 
it likely that a range of cut-off values exist. In an equi- 
librium situation synchrotron and inverse Compton losses 
will steepen the spectrum by a Aq — —1 (Longair 1994), 
but as soon as the electrons are sufficiently far removed 
from the shock front the spectrum is more likely to have 
an exponential cutoff. An exponential cut off is also more 
likely, if the spectrum near the shock front is age limited 
instead of loss limited (Reynolds 1998). 

The rims have a typical width of 1.5"-4" and clearly 
stand out above the background and diffuse emission 



(Fig. 5). We can estimate a typical loss time from the 
widths by combining them with the plasma velocity away 
from the shock front. The shock velocity in Cas A has 
been estimated from proper motion studies using Ein- 
stein and ROSAT HRI X-ray images, and is thought to 
be Vs = 5200 km s^^ (Vink et al. 1998). The velocity of 
plasma away from the shock front, for a standard shock 
compression ratio of 4, is 0.25us = 1300 km s~^. This 
value is consistent with X-ray Doppler shifts, which imply 
a velocity relative to the shock front of Av ^ 1400 km/s 
(Willingale et al. 2001a). For a distance of 3.4 kpc the 
measured widths of the rims correspond to a typical loss 
time of 18 to 50 yr. 

The energy history of an electron with initial energy i?o, 
suffering synchrotron losses, is: 

^^'^ = l + Eom/635- 
After a time r (eq.l), the electron will have an energy ^Eq, 
whereas higher energy electrons, after the same amount of 
time, will have energies E'{t) < Eq. The relation between 
the peak emitting frequency and electron energy is given 
by (Ginzburg & Syrovatski 1965): 

Vpeak = 1.8 X IQ'^B^E''. (4) 
In reality we have to take into account a range of frequen- 
cies, as the the synchrotron emission function is broad. 

The rim is best visible in the X-ray continuum domi- 
nated band of 4-6 keV, and in some narrow bands in be- 
tween line emission features. Taking 5 keV (1.2 x 10^^ Hz) 
as the typical photon energy, we can combine eq. 4 with 
eq. 1 to obtain a typical electron energy of = 65 erg 
to 91 erg and B = 0.16 mG to 0.08 mG. The power slope 
of the Chandra spectrum near the shock is —2.2, which 
indicates that the photons are emitted by electrons close 
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to the break in the electron distribution. In fact using 
the the srcut-model (cf. §3), we derive for B = 0.1 mG, 
Ebreak ^ 25 erg. However, the spectrum is hkely to be 
an integration of spectra with different break energies. So 
close to the shock front Ebreak is likely to be higher. 




5 -5 

Angular distance (arcsec) 

Fig. 5. Emission profiles of the rim in the southeast (left) 
and northeast. The profiles were extracted using strips of 20 pix- 
els wide (9.8") from a Chandra image, which was a combination of 
continuum dominated bands between 2.7 and 9 keV. 



In Fig. 6 we have illustrated our estimates of magnetic 
field strength and electron energy graphically. Apart from 
synchrotron losses we also have included inverse Comp- 
ton losses. For this we have estimated the photon energy 
density from the cosmic microwave background and the 
far infrared emission from Gas A itself, the latter being 
the most dominant factor. From the IRAS infrared fluxes 
listed by Braun (1987) we estimate a photon energy den- 
sity of M = 4.7 X 10^^^ ergcm^'' , corresponding to an 
equivalent magnetic field of i? = 1.1 x 10~^ G. 

The estimates of the magnetic field near the shock front 
give values of roughly a fifth of the average magnetic 
field strength and indicate a preshock magnetic field of 
~ 2 — 4 X 10~^, a factor ten higher than the canonical in- 
terstellar medium value. This may cither hint at a higher 
magnetic field surrounding Gas A, for instance because 
the wind of the Wolf-Rayet progenitor carried a high mag- 
netic field, or the presence of accelerated particles near the 
shock front has given rise to magnetic field amplification 
due to non-linear plasma wave growth. The first hypoth- 
esis was suggested by Biermann & Gassinelli (1993) in or- 
der to explain the acceleration of cosmic ray by supernova 
remnants up to 10^^ eV for which high magnetic fields are 
necessary. The second hypothesis was suggested by Lucek 
& Bell (2000) for identical reasons. 
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Flc;. 6. — The shaded area indicates the most likely values of Smaa; 
and B just behind the shock front. This area is formed by the two 
dashed lines, which mark the electron energies that contribute to the 
continuum emission around 5 kcV, allowing for emission at l/3rd 
and 10 times the peak frequency for synchrotron emission, and the 
solid lines, which arc formed by assuming the width of the rims is 
determined by radiative losses (decay times of 18 to 50 yr and a 
shock velocity of 5200 km/s arc assumed). The width of the rims 
should at least be larger than twice the gyro-radius. This means 
that for a typical width of 3", the hatched region is excluded. 



6. THE MAGNETIC FIELD STRENGTH AND THE NATURE 

OF THE TEV EMISSION 

The measured magnetic field strengths have important 
consequences for the gamma ray emission from Gas A. 
For instance, our lower limit on the average magnetic 
field strength implies that the gamma ray emission from 
bremsstralilung is almost an order of magnitude below the 
CGRO-EGRET upper limit. 

It is also understood that for the high magnetic field 
in Gas A, it is unlikely that the detected TeV emission is 
inverse Gompton emission (Aharonian et al. 2001; Atoyan 
et al. 2000). We have investigated this important issue by 
constructing a two zone model. Zone 1 is responsible for 
the bulk of the radio emission, such as the bright radio 
ring. It should have B > 0.5 mG according to lower limit 
on the magnetic field. As the electron population in this 
zone is old, no X-ray synchrotron emission is produced. 
We adopt a somewhat arbitrary maximum electron energy 
of 0.1 TeV, but in reality a range of values exist inside the 
remnant The inverse Gompton emission from this zone is 
negligible. 

Zone 2 is responsible for the X-ray synchrotron emis- 
sion from the rim and has B ^ 0.08 — 0.16 mG and 
Emax ~ 40—60 TeV. Based on our flux measurements, and 
assuming spherical symmetry, we estimate that the X-ray 
synchrotron emission should be less than one third of the 
X-ray emission around 5 keV, but as we focused on a par- 
ticular bright part of the rim, 1/lOth to l/20th seem more 
likely values. Moreover, for B = 0.08 mG and an X-ray 
synchrotron flux contribution of l/3rd, the non-thermal 
bremsstrahlung flux at 100 keV is above the measured up- 
per limit. This criterion is relaxed for higher magnetic 
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fields near the shock front. 
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Fig. 7. — A simplified two zone model of the broad band radiation 
from the cosmic ray electron population, i.e. synchrotron radiation 
(dotted line), brcmsstrahlung (dashcd-dotted), and inverse Comp- 
ton emission( dashed line). The solid line gives the total emission 
from both zones. The parameters near the shock front (zone 2) 
were chosen to be B = 1.0 x IQ-* G, Emax = 55 TeV. B = 0.5 mG 
for zone 1. The data points give the flux at 1 GHz, measurements 
by ?, X-ray continuum emission {ASCA and BeppoSAX), CGRO- 
COMPTEL (Strong et al. 2000), the CGRO-EGRET upper limit 
(Esposito et al. 1996), and HEGRA (Aharonian et al. 2001). 

This two zone model is different from the two zone model 
proposed by Atoyan ct al. (2000), for which the two zones 
relate to the radio knots and diffuse radio emission. We 
adopt here the view of Anderson & Rudnick (1996), that 
there are no intrinsic differences between the electron cos- 
mic ray populations of the radio knots and the diffuse 
emission, except that in the radio knots the magnetic field 
is stronger. Note that the electron cosmic ray density is 
higher in zone 2. Zone 1 is brighter at radio wavelength 
due to the stronger magnetic field. 

In order to estimate the inverse Compton emission we 
need to estimate the photon density caused by Cas A's far 
infrared emission. IRAS flux measurements (Braun 1987) 
indicate a total photon flux of Fp/j = 3.9 x 10^ ph s~^ 
cm~^, centered around 10^^ Hz. The far infrared emission 
is associated with ejecta (Lagage et al. 1996). Approxi- 
mating the geometry by a sphere with radius similar to 
the ejecta, r = 1.8 pc, and assuming that half the number 
of photons spend an average time of r/c inside the sphere, 
we obtain for the photon density^: 

n,, = ^^^ = 72cm-3, (5) 

where d = 3.4 kpc is the distance to Cas A. A similar rea- 
soning gives a somewhat lower photon density for zone 2, 
but the photon field is much more anisotropic, on average 
coming from the direction of the remnant's center. For the 
spherical symmetry assumed, the anisotropy has no influ- 
ence on the total inverse Compton emission, but as the 
maximum upscattering occurs for head on collisions, most 
energetic inverse Compton emission from zone 2 will be 
directed toward the center of the remnant, giving rise to a 
centrally brightened morphology, despite the fact that ac- 
cording to our model the dominant inverse Compton emis- 
sion is coming from a narrow shell near the shock front. 

^ The other half of the photons escape directly, but all photons will 



In contrast, brcmsstrahlung from the same electron pop- 
ulation will show a ring-like morphology. Unfortunately, 
these interesting differences in morphology occur on the 
scale of arcminutes, which is currently beyond the reach 
of gamma-ray instruments. 

Our two zone model is illustrated in Fig. 7, employ- 
ing typical values for B and Ebreak near the shock front. 
The inverse Compton emission is completely dominated by 
the electrons in zone 2. In Fig. 7 inverse Compton emis- 
sion is well below the TeV flux by HEGRA. It is there- 
fore likely that the TeV emission is dominated by pion 
decay. However, we cannot yet make a definitive state- 
ment regarding the inverse Compton emission. For lower 
values of the magnetic field zone 2 is also the dominant 
source of brcmsstrahlung by relativistic electrons. If the 
magnetic field is at the low end of our inferred magnetic 
field, B = 0.08 mG, and the far infrared photon density 
is larger than 150 cm^'', the HEGRA detection, and the 
BeppoSAX/ GGRO-OSSE upper limit on the non-thermal 
bremsstrahlung arc still consistent with inverse Compton 
emission. The basic uncertainty in the photon density 
stems from the uncertain geometry, which for zone 2 may 
also include anisotropic photon emission. The combina- 
tion of a low magnetic field and a rather high infrared 
photon density seems unlikely, but is not impossible. Even 
lower magnetic fields in zone 2 are unlikely, as the non- 
thermal brcmsstrahlung would exceed the upper limit re- 
ported in §3. 

7. CONCLUSION 

We have provided estimates of the average magnetic 
field and the magnetic field near the shock front. The 
lower limit to the average magnetic field is based on the 
assumption that the electron spectrum follows a power law 
distribution in momentum, or is steeper than that at en- 
ergies between 100-1000 keV. 

The magnetic field strength near the shock front, was 
measured from the width of the narrow rim surrounding 
Cas A in X-ray continuum images, assuming that the emis- 
sion mechanism is synchrotron emission, and the width is 
determined by synchrotron losses. The contribution of the 
X-ray synchrotron emission to the X-ray continuum flux 
is less than one third. This is supported by the X-ray 
morphology, which shows that even around 10 keV a sub- 
stantial part of the continuum is associated with the bright 
ejecta shell (Bleeker et al. 2001). 

This justifies taking into consideration alternative mod- 
els than synchrotron radiation for the hard X-ray emis- 
sion of Cas A. This is was done in §2, where we tested the 
lower hybrid wave model. The brcmsstrahlung spectrum 
predicted for electrons accelerated by lower hybrid waves 
has a similar spectral slope as is observed. The maxi- 
mum energy to which the electrons are accelerated by this 
mechanism is ~95 keV, according to fits to the BeppoSAX 
data. 

The estimated lower limit on the average magnetic field 
is 0.5 mG, whereas we find ~ 0.1 mG near the shock front. 
So even near the shock front the magnetic field strength is 
an order of magnitude above what is expected from shock 
compression of the canonical interstellar magnetic field of 
3 /xG. This can be explained by either a higher magnetic 

to pass through zone 2. 
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field strength in tlie medium surrounding Cas A, possibly 
as the result of the stellar wind of the progenitor (Bier- 
mann & Cassinelli 1993), or it means that the magnetic 
field is being rapidly enhanced near the shock front (Lucek 
& Bell 2000). 

Electron energies close to the shock front may reach 
~50 TeV or higher, but in contrast to Reynolds & Keohane 
(1999), we suggest that most of the cosmic ray electrons 
responsible for the bright radio ring are in an environment 
with a high magnetic field, where due to synchrotron losses 
the the cut-off energy is well below the energy where X-ray 
synchrotron emission can be expected. 

The high average magnetic field makes it likely that TeV 
emission from Cas A is dominated by pion decay (Atoyan 
et al. 2000), but there is some uncertainty, as in our two 
zone model of Cas A it is still possible to obtain a high 
inverse Compton flux by assuming that the magnetic field 
is at the low end of what we have inferred from the X-ray 
rims, B = 0.08 mC, and using a higher far infrared photon 
density of 150 cm~^. 

Future mission will improve on those results. For 
instance, GLAST will probably be able to detect 
bremsstrahlung in the GeV range, which will provide a 
much more direct way of estimating the average magnetic 
field of Cas A. 

As for the hard X-ray emission. INTEGRAL will not 
greatly improve the current measurement of the hard X- 
ray continuum emission, but by accurately measuring the 
'*'*Ti related line emission, it will constrain that component 



of the emission around 80 keV, and thus indirectly im- 
prove our understanding of the continuum emission around 
100 keV. 

One of the consequences of our analysis is that the 
gamma-ray emission will be mostly associated with the 
outer rim of Cas A, not with the bright X-ray and ra- 
dio shell. Unfortunately, the spatial resolution required to 
measure this will not be available in the near future. 

Finally, as an interesting by product of our analysis, we 
mention that the new modeling of the hard X-ray emis- 
sion with the lower hybrid model yields revised values of 
the ^^Sc line emission at 68 keV and 78 keV, which cor- 
respond to a 50% higher mass of ^^Ti synthesized in the 
explosion compared to what was reported by Vink ct al. 
(2001), i.e. M(44Ti) = 1.8 x 10"* M©. This is closer to 
the latest CGRO-COMPTEL measurements (Schonfelder 
et al. 2000). 
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edge the use of the HEASARC (Goddard Space Flight 
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APPENDIX 

APPENDIX 

In this appendix we want to demonstrate how free parameters in the fit of the model continuum (the maximum 
accelerated electron energy and the fraction of ambient electrons that are accelerated) are in principle determined by 
the plasma physics parameters of Cas A. The maximum electron energy is related to the maximum electric field that a 
lower hybrid wave can develop and the time that an electron can spend in this field, governed by the thickness of the 
shock precursor. The fraction of electrons that may be accelerated is expected to be determined by a marginal stability 
criterion, i.e. the electrons will be accelerated until the heating of the ambient plasma produced by Coulomb colisions of 
the accelerated electrons violates the conditions necessary for lower hybrid waves to exist, i.e. that the electron gyroradius 
be less than the wavelength divided by 27r. 

The lower hybrid waves are generated by cross field ion motions produced by shocks. A fraction of preshock ions are 
reflected back upstream ahead of the shock with a velocity relative to the preshock plasma of ~ 2vs, where Vs is the 
shock velocity. Ions reflected initially back along the shock velocity vector at a perpendicular shock travel a distance 
0.68us/r2i before returning (Gedalin 1996). Here f7j is the ion gyrofrequcncy. The maximum distance ahead of the shock 
that reflected ions can reach \s d = LZlVs/^i- This obtains for ions reflected from the shock at an angle 30° to the shock 
front, so that 2vs sin 30° = Vs- 

Consider a shock front at an angle 9 to the local magnetic field direction. In order for the lower-hybrid waves generated 
ahead of the shock by the reflected ion precursor to grow to large amplitudes we require (Laming 2001b, i.e. the wave 
group velocity away from the shock is equal to the shock velocity, so that the time available for wave growth is essentially 
unlimited on the timescale of these processes) 

du _ Vs _ auj , , 

dki. ~ cose ~ k± ^ ' 

where a = 1/ (2 cos /3 cos 0) and /3 is the angle between doj/dk and the reflected ion bulk velocity vector U. A minimum 
requirement for wave growth is —1 < a < —1/2. Since it is ions returning to the shock that can generate the necessary 
waves, waves generated at the earliest possible period in the reflected ion orbit, i.e. lowest value of |cos/J|, when a — 1, 
(for cos^ ~ 1) will dominate the growth. 

The maximum electric field that the lower hybrid wave may develop before ion trapping and heating set is given by 
Karney (1978, in mks units); 



' w y 4A;_L V w / 4cos6' ^ ' 
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Taking the wave frequency lu — ^lh = v^W^, and the component of the electric field directed along the magnetic field 
direction to be Eupi/iVpe gives a ~ UlnVs/^QcosO for electron acceleration along the magnetic field direction. Qlh is 
the geometric mean of the ion and electron gyrofrequencies and is known as the lower hybrid frequency. 

We treat the case of a nonrelativistic shock which may accelerate electrons to relativistic energies. We work in the rest 
frame of the prcshock medium. In this frame an electron is accelerated by the component of the electric field directed 
along the magnetic field, until it outruns the reflected ion shock precursor a distance d from the shock front. Assuming 
acceleration begins with the electron at the shock front moving along B with initial velocity vq, it will continue for a time 
t given by 



a \ 



(at + vo/\/l-v^/cA 2 
l + - -n 



> sin0 = d + Vst. (A3) 



which gives a quadratic equation for t: 

(c^ sin^ 6-v',)+tl (^;o sin^ 6 - sin 6) - 2dv, ] - - ^ff''"f^^ = (A4) 

We require sm9 > Vg/vo in order that the electron stays ahead of the shock front, and maximum acceleration time will 
be realized for sin^ = Vs/vq. In this case cos^ ~ 1 and the term in t in equation A4 simplifies to —2dvst and 



d 

t = 



l + (cV^^-l)ll + - ?^L==| S. (A5) 



For Vq « c, at = ^/2advoJv^ and for vq c, at = 2ad/vs (c^/uq — l). Since ad a it is easy to see that the highest 
electron energies are most efficiently reached by fast shocks in hot plasma, i.e. for the highest plausible values of Vg 
and Vq in Cas A. At sin0 = Vs/vq, three to four shock interactions are required to accelerate electrons to the observed 
maximum energy of ~ 95 keV, increasing to about ten shock interactions if sin 61 = 1.2vs/vo, for Vg = 3000 km s""'^ and 
Vq = 3x 10^ cm s~^ (corresponding to an electron temperature of 4 x 10^ K). Even faster shock require fewer interactions. 
The maximum energy in the accelerated electron distribution will be determined by the number and obliquity of shock 
interactions experienced by these electrons during the evolution of Cas A, and so for the time being remains uncertain. 
However the maximum energy inferred from Fig. 1 is certainly plausible in these terms. In principle, given enough shock 
interactions, relativistic electrons can be produced. If a means of transporting these electrons to the forward shock were 
available, this might have a bearing on providing an injection mechanism for further electron acceleration by a Fermi 
mechanism. 

The normalization of the non-thermal bremsstrahhmg spectrum is determined by the fraction of plasma electrons that 
are accelerated in this way, and this will be determined ultimately by a marginal stability criterion. The instability 
requires the electron gyroradius to be less than the wavelength divided by 27r, and once a sufficient number of electrons 
are accelerated to heat the ambient plasma to temperatures higher than this, the instability will shut down. The ambient 
plasma loses energy through radiation, adiabatic expansion and possibly conduction. In Laming (2001b) the heating rate 
due to the observed accelerated electrons is estimated to be 4 x IQ^^^rieTi'^ ergs cm~"^s~^ for ambient and accelerated 
electron densities Ug and n'^. The electron heating rate is dT/dt = 0.019ng K s~^, and equating this to the cooling rate by 
adiabatic expansion of the Cas A remnant yields an accelerated electron fraction rig/rie ~ 0.02. The inclusion of ionization 
and radiative power losses and relaxing the steady-state assumption inherent in this estimate will allow higher accelerated 
electron fractions. 

In Laming (2001b) an accelerated electron fraction of 4% was found to give a good match to the ~ 50 ksec of BeppoSAX 
MECS and PDS data on Cas A then available. The estimate herein using only the ^ 500 ksec PDS data is higher at 
around 10% (see Fig. 1), but a range of plasma temperatures exists in Cas A (c.f. Vink et al. 1996; Willingale et al. 2001a), 
and only the hottest plasma with a temperature of ^ 4 x 10^ K contributes to the PDS bandpass. So with respect to all 
the hot plasma, an accelerated electron fraction of ~4% is consistent with the data, and preferred on theoretical grounds. 
Simulations of the evolution of the electron temperature following electron heating indicate that such an accelerated 
electron density heats plasma in the location of the contact discontinuity to temperatures in the range 4 — 5 x lO'' K. 
These temperatures are basically constant over the past ~ 30 years over which X-ray observations have been made. In 
this respect simulations of Cas A expanding into an p cx 1 /r^ circunistcilar medium appear to represent the physics much 
better than the simulations using a uniform density circumstellar medium in Laming (2001b). 
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Table Al 

Summary of model fits to the BeppoSAX-P'DS spectrum. 





LHW model 


srcut 


sresc 


power law 


Norm continuum'^ 


0.058 ± 0.002 


2522 


2522 


2.37 ±0.2 


power law index 








-3.32 ± 0.05 


Emax/Scuto// (keV) 


95 ± 5 


1.21 ±0.01 


3.42 ± 0.02 




Norm thermal emission^ 


0.62 ± 0.02 


7.6 ±3.3 


0.0(< 2.5) 


0.0(< 17) 


"Sc line flux (lO-*^ phcm-^s-^) 


3.2 ± 0.3 


2.7 ± 0.4 


2.9 ± 0.4 


2.1 ±0.4 


Ta K(v flux (10"'' phcm^-s^M 


2.4 J- 0.7 


2.6 + 0.8 


2.8 + 0.8 


1.4 + 0. 




27:i,l/2ol 


27:!.()/2.".2 


28y,:i/252 


2(i7,:;/251 



^ The normalization for the LWH model is defined as 10^^ nuneV /47rd'^ cm~^. The normaliza- 
tion for the synchrotron models is fixed to the radio flux density at 1 GHz for the epoch 2000.0, 
and radio spectral index is 0.775. These values are based on a compilation of radio data, and 
includes the effect of the secular flux decrease of Gas A (Baars et al. 1977). The normalization for 
the power law model is ph/s/keV at 1 keV. 

The normalization is defined as 10^^ nun^-V / And^ . The temperature was fixed to 4 X lO"^ K. 
Note that this definition assumes solar abundances, which is not applicable for Gas A, so care 
should be taken before using it to calculate electron densities. 
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Table A2 

Best fit parameters for spectra from the northeastern rim. 



NE Rim Inside NE Rim 



kT (keV) 


3.7 ±0.7 3.0 ±0.4 


net (10^° cm-^ s-^) 


3.3 ±0.5'' 5.2 ±0.6 


Ne 


4+^ ° 


Mg 


0.9 ± 0.4 


Si 


4.5 ± 0.6 


S 


4.3 ± 0.7 


Ar 


4.4 ± 1.5 


Ca 


5.4 ±2.5 


Fe 


2.2 ±0.4 


EM ■=(10^° cm-S) 


3.0 ±0.4 7.5 ±1.0 


PL index 


-2.24 ±0.05 -4.5(< -4.3)"^ 


PL norm (10"-'' s^^cm"^ ") 


2.5 ±0.2 3.4 ±0.6 


PL contribution 4-6 keV 


84% 11% 


Nh (10^^ cm"^) 


1.08 ± 0.03 




261/161 363/153 



Note. — The spectral model consisted of SPEX^s non- 
equilibrium ionization model (vl.lO), and a power law (PL) con- 
tinuum, in order to approximate the synchrotron radiation. The 
column density and abundances were first determined for the spec- 
trum inside the NE Rim, and subsequently applied to the rim spec- 
trum. Following Vink et al. (1996) wc set He and N abundances to 
10 times solar (Anders & Grevesse 1989). The poorly constrained 
O abundance was coupled to Si; Ni was coupled to Fe. Errors 
correspond to 95% confidence intervals. 

^Note that this value is roughly consistent with the synchrotron 
loss times, if the electron density is Ue ~ 30 cm~^. 

^The He and N abundances were assumed to be 10 times solar, 
following Vink et al. (1996). 

^The emission measure (EM) is defined as nurieV /'^Trd^ . 

'^The power law slope was poorly constrained, but tended to 
be steep, we therefore fixed it to a value of 4.5. The lower limit 
corresponds to Ax^ = 4. 



^At 1 keV. 



